www.gusucode.com > 信道估计算法源码程序 > 信道估计算法源码程序/CS-Channel Estimation/MMSE_MSE_calc.m

    %Function Declaration:
function ms_error=MMSE_MSE_calc(X,H,Y,Rgg,variance,N,L);
%This function generates mean squared error for the the MMSE estimator..
%EVALUATION OF Hmmse
%Hmmse=F*Rgg*inv(Rgy)*Y;
F=fft(eye(N,L));
I=eye(N,N);
Rgy=Rgg * F'* X';
Ryy=X * F * Rgg * F' *X' + variance * I;
for i=1:N
    yy(i,i)=Y(i);
end
Gmmse=Rgy * inv(Ryy)* Y;
Hmmse=F*Gmmse;

ms_error_mat=mean(((abs(H)-abs(Hmmse))/abs(H)).^2);
for i=1:N
    if(ms_error_mat(i)~=0)
        ms_error=ms_error_mat(i);
    end
end